Drag force on an oscillating object in quantum turbulence 
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This paper reports results of the computation of the drag force exerted on an oscillating object in quantum 
turbulence in superfluid ^He. The drag force is calculated on the basis of numerical simulations of quantum 
turbulent flow about the object. The drag force is proportional to the square of the magnitude of the oscillation 
velocity, which is similar to that in classical turbulence at high Reynolds number. The drag coefficient is also 
calculated, and its value is found to be of the same order as that observed in previous experiments. The corre- 
spondence between quantum and classical turbulences is further clarified by examining the turbulence created 
by oscillating objects. 

PACS numbers: 67.25.dk, 47.37.+q 



I. INTRODUCTION 

In recent years, quantum turbulence has become an estab- 
lished research area in low-temperature physics.''^ The sim- 
ilarities and differences of quantum turbulence in relation to 
classical turbulence comprise one area of interest of quantum 
turbulence. For example, although both classical and quantum 
turbulences consist of eddies or vortices, these two kinds of 
turbulences are quite different in terms of their hydrodynamic 
properties. 

First, classical fluid can have arbitrary vorticity, while in 
contrast the vorticity in superfluid ^He is quantized; circula- 
tion around a quantized vortex isK — h /m, where h is Planck's 
constant and m is mass of ^He atom, with the core size of 
a quantized vortex given by the healing length of superfluid 
^He, which is about the atomic size of ^He. Moreover, quan- 
tized vortices cannot be generated under usual experimental 
conditions, unlike eddies in classical fluids. Instead, remnant 
vortices are essential for the generation of quantum turbu- 
lence. It is believed that remnant vortices are usually pinned to 
the roughness of the wall.^ Indeed, Hashimoto et al^ experi- 
mentally showed that quantum turbulence cannot be generated 
without remnant vortices. 

Second, classical fluid is viscous, while superfluid ^He con- 
sists of viscous normal fluid and an inviscid superfluid ac- 
cording to a two-fluid model. The two fluids are coupled in 
turbulence through mutual friction, which is the interaction 
between vortices and normal fluid. The ratio of the two flu- 
ids is temperature dependent, and the normal fluid dominates 
the fluid if the temperature is close to the transition temper- 
ature, while in the limit of zero temperature, the superfluid 
component dominates the fluid, and pure quantum turbulence 
is achieved. The attention of the present study is confined to 
this latter limiting case. 

In spite of the above distinct properties, quantum and clas- 
sical turbulences also exhibit similarities. Kolmogorov's law, 
which is a statistical law found in a classical turbulence, has 
also been observed in quantum turbulence both experimen- 
tally and numerically i^i^i^i^i^ This law can thus be taken to be 
universal over quantum and classical turbulences. 

Under this new motivation, quantum turbulence has re- 
cently been studied, and accordingly the progress of exper- 
imental techniques has given rise to a variety of methods 



to create and observe quantum turbulence. To both create 
and detect quantum turbulence, objects oscillating in a super- 
fluid, such as spheres, grids, thin wires, and tuning forks are 
often employed. The most commonly mea- 
sured physical quantity is the velocity of the oscillating ob- 
jects. At temperatures lower than 1 K, the velocity exhibits 
the clear transition from laminar to turbulent state when the 
driving force is increased.— For magnitudes of velocities less 
than about 50 mm/s, the drag force is proportional to the os- 
cillation velocity due to collisions between the oscillating ob- 
ject and some remaining excitations. At higher magnitudes of 
velocity, the drag force is proportional to the square of the os- 
cillation velocity magnitude because turbulence is generated 
around the object. 

These results also describe a similarity between quantum 
and classical turbulence; in classical fluid, the drag force on an 
object in a uniform flow at high Reynolds number is described 
by 



(1) 



where Co is the drag coefficient, p is the density of the fluid, 
A is the projection area of the object normal to the flow, and 
U is the flow velocity.'^ This relation could be applicable 
to quantum turbulence created by an oscillating object. At 
low Reynolds number, Stokes's drag force acts on the object, 
which is proportional to the magnitude of the flow velocity, 
with the result that the drag coefficient Co becomes inversely 
proportional to the magnitude of the flow velocity. Fitting ex- 
perimental results to a relationship of the form of Eq. ([T]), we 
derive a constant drag coefficient, which is of order unity, as 
for classical turbulence. 

In this paper, we shall further add to knowledge of the sim- 
ilarities of quantum and classical turbulences by presenting 
numerical simulations of the drag force and coefficient on the 
oscillating object. However, in calculating the drag force nu- 
merically, the pressure on the object needs to be calculated 
over the entire surface of the object. The pressure at each 
point on the object can be calculated from the velocity field 
by Euler's equation, and the velocity field is generated by the 
vortices. It is possible to calculate the pressure if the num- 
ber of the vortices is relatively small, but it is difficult with 
current computational resources to integrate the pressure in 
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quantum turbulence in cases where a large number of vortices 
exist. As a consequence, we shall introduce a method to eval- 
uate the drag force. We assume that the oscillating object does 
work on the fluid when the vortices are grown by the object 
and that a drag force is caused as a result of the work. It is 
shown later that the drag force can be estimated from the en- 
ergy given to the fluid by the object in an equilibrium state 
of quantum turbulence. The increase in the kinetic energy of 
the fluid can be evaluated from the energy of vortices grown 
by the object. Finally, from the energy given per unit time 
divided by the velocity of the object we derive the drag force. 

In this paper, the drag force in quantum turbulence is tar- 
geted for the following reason: the drag force becomes defi- 
nite according to the method of this paper only if the vortices 
are grown by the oscillating object. In the laminar regime, it 
is expected that the drag force is caused by collisions between 
the object and excitations and that the growth of vortices is not 
essential. In the turbulent regime, on the other hand, the drag 
force is caused by the growth of vortices. For this reason, we 
focus on the dynamics of vortices around an oscillating object 
and solve the equations of motion for the vortices numerically. 

In the simulation, we consider an object oscillating with a 
constant magnitude of velocity v, and estimate the drag force 
due to the vortices, that is the drag force Fdrag is described 
in the form F^rag = f{v)- Note that this is in contrast to ex- 
perimental studies, which measure the oscillation velocity of 
the object for a certain driving force, obtaining the velocity 
V — (-Pdrag)- The functional relation Fdrag = f{v) used in 
the present simulations is thus inverse to the relation used in 
the experiments, but it is worth comparing the results for the 
following reason. As we shall see in Sec. Hill when the ob- 
ject oscillates, remnant vortices are stretched and the vortex 
tangle grows, which results in the drag force on the object. Fi- 
nally, the vortex tangle reaches an equilibrium state of quan- 
tum turbulence. The drag force is then uniquely determined 
by the oscillation velocity in the simulation. In an experimen- 
tal equilibrium stat e on the other hand, the drag 
force responsible for the dissipation by vortices must be equal 
to the driving force that injects the energy into the system. 
In addition, the oscillation velocity is measured under a con- 
stant driving force, namely, drag force in an equilibrium state 
of turbulence in the experiments. Hence there is one-to-one 
correspondence between the oscillation velocity and the drag 
force, which allows us to compare the drag force of the simu- 
lation with that of the experiments. 

The coupled dynamics of vortices and the oscillating object 
is not considered in this work, which is validated by verifying 
that the drag force depends only on the oscillation velocity. 
The oscillation of the object causes the growth of vortices, 
and the rate of the growth depends on the oscillation velocity 
because the vortices attached to the object are stretched by the 
object (details in SecHIIll. The vortex dynamics are influenced 
only by the magnitude of the oscillation velocity irrespective 
of whether the dynamics of the object are considered or not. 
This is because since the drag force is evaluated from the rate 
of the vortex growth, the drag force does not depend on the 
details of the motion of the object but the oscillation velocity. 

A previous numerical simulation of an oscillating sphere 



was performed by Hanninen et al. They used the same size 
of the sphere and frequency of the oscillation as the exper- 
iment of Shoepe et aly— The simulation seems to show a 
growth of vortices toward turbulence, but an equilibrium state 
of quantum turbulence is not achieved, which differs from the 
present results. This is because in Hanninen et aZ.'s'^ study 
the vortices extending between the wall and the sphere re- 
mained attached and continued to generate vortex rings, while 
in our simulations, mature vortices are soon detached from 
the object and the remaining vortices are able to grow succes- 
sively (details in SecHIlli. It is important to note that the drag 
force cannot be calculated from the simulations of Hanninen 
et al}'^ in the manner of the present study, since our method 
assumes an equilibrium state. 

The contents of this paper are as follows. In Sec. [Ill we 
shall clarify the configuration and formulation of the model 
and introduce the equations of motion required in the dynam- 
ics of quantum turbulence. In Sec. Hill we shall describe the 
process whereby quantum turbulence is generated by an os- 
cillating object. In Sec. Hyl we shall introduce the present 
proposed method of evaluating the drag force on the oscillat- 
ing object and show the dependency of the drag force on the 
velocity of the object. Section |Vl is devoted to the conclusion 
and the discussion. 



II. MODEL AND FORMULATION 

Superfluid ^He at K can be treated mathematically as an 
ideal incompressible fluid of vanishing viscosity. In a super- 
fluid, any circulation is quantized in units of K. The vorticity 
is infinitely distributed on thin vortex core, whose size is the 
order of 1 A, so that the circulation around the vortex filament 
is restricted to K.^ With these features, it is valid to adopt the 
vortex filament approximation, which regards a vortex as a 
line. Numerical study of quantum turbulence can be accom- 
plished by solving the equations of motion for the vortex fila- 
ments that form the turbulence. According to Helmholtz's 
theorem, a vortex filament moves with superfluid velocity on 
the vortex. In an infinite system of superfluid, the superfluid 
velocity is equal to the velocity generated by vortex filaments, 
which takes the form of the Biot-Savart integration, 

'"u>{r)^— — ■ ^ — , (2) 

where Si is a position vector on a vortex and the integration is 
performed over all the vortex lines. In attempting to determine 
a velocity field at a certain point r = s on a vortex, Eq. ^ 
diverges as Si s. The divergence is evaded by introducing 
cutoff parameter aq in Eq. (|2]i. 

In the presence of a boundary, additional superflow Vb ap- 
pears so as to satisfy the boundary condition for an inviscid 
fluid, 

• n = 0, 

Vs=v^ + vi„ (3) 

where n is a vector normal to the boundary. While many pos- 
sible geometries may be considered for the vibrating objects. 
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FIG. 1 : Image vortices in a sphere for a vortex element. 



in the present study we restrict attention to the sphere for sim- 
plicity. 

For a velocity field generated by an infinitesimal element 
ds of a closed-loop vortex, Eq. Q is satisfied by considering 
an image vortex inside the sphere. ^^'^^ The image vortex con- 
sists of two parts: a radial image ds' r and a tangential image 
ds't. In Fig.[Tl the origin is set at the center of the sphere, a is 
the radius of the sphere, s is a radius vector of the vortex ele- 
ment, dsr is a radial component of ds, and dst is a tangential 
component defined as dst ^ ds — dSr- The tangential image 
is located at a^/s from the origin and anti -parallel to dst- Its 
length is {a/s)^dst, and it has circulation {s/a)K. The radial 
image directs from ds'f. to the origin, its length is a'^/s and it 
has circulation {dsra)K. 

When the sphere moves with the velocity Up, the boundary 
condition [Eq. Q] is satisfied in a reference frame that moves 
with the sphere: 

{vs - tip) • n = 0. (4) 

The extra term in the parenthesis is canceled out by adding the 
following term to Vf,: 

Vuir) = V<I>„(r), 

where is a flow induced by the motion of the sphere and de- 
termined by the instantaneous velocity of the sphere. Thus 
Vh is obtained by adding the contribution from the image vor- 
tices and the velocity field given by Eq. (|5]l. 

The dynamics of the vortices are completely determined by 
the above formulation. When two vortices cross or when a 
vortex becomes close to a spherical boundary, it is theoreti- 
cally known that reconnection, which is topological change, 
occurs. In the framework of the vortex filament model, the 
details of the vortex core are neglected, and so these formu- 
lations do not describe the dynamics of the reconnection. In 
the present numerical simulation, an exceptional routine is in- 
cluded in our numerical code to allow reconnection to occur 
when a vortex approaches another vortex or a spherical bound- 
ary within the numerical resolution. 

It is thought that quantum turbulence originates from rem- 
nant vortices rather than from the intrinsic nucleation of 
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FIG. 2: (Color online) The time evolution of turbulence generation 
for the case of a sphere oscillating with a velocity magnitude of 90 
mm/s. See the text for further details. 

vortices. ^-^ '^ Thus we inject vortex rings toward the oscillat- 
ing object to induce quantum turbulence, considering the ex- 
periment of the Osaka group^i^ This approach using the injec- 
tion of vortices requires the consideration of two parameters: 
the ring size and the time interval in which the vortices are 
injected. The ring size can be estimated experimentally by the 
time of flight for the vortex ring, and we utilize the value of the 
ring size 1 /jm.'^ We assume that the time interval is smaller 
than the period of the oscillation, say, 0.05 ms.'^ In order to 
replicate the conditions of Ref. 15, in the present study the 
same parameter values as this previous study are used, such 
as the radius of sphere and the oscillation frequency of the 
sphere. The diameter of the sphere is 3 /jm, the frequency 
of the oscillation is 1590 Hz, while the oscillation velocity is 
chosen in the range of 30-90 mm/s. 

III. DEVELOPMENT OF TURBULENCE 

The time evolution of the turbulence generation is shown 
in Fig. |2] and Ref. 27. The sphere oscillates horizontally, and 
vortex rings are injected from the bottom of the medium [Fig. 
|2ja)]. When the vortex rings collide with the sphere, recon- 
nection occurs and the vortices are attached to the sphere. 
It can be seen that the attached vortices are stretched as 
the sphere moves [Figs. |2jb) and|2|c)]. Due to the succes- 
sive injection of vortex rings, the process is repeated and 
the stretched vortices form a tangle around the sphere [Fig. 
|2jd)]. The vortices grow in size and are then detached from 
the sphere as follows. The flow caused by the motion of the 
sphere [Eq. ©J drives the end points of the attached vortices 
to the stagnation point of the sphere. A pair of the end points 
of the vortex approaches each other as they converge to the 
stagnation point. Finally, reconnection of the pair of the end 
points then occurs and the vortex is detached from the sphere 
[Fig- Etc)]. In spite of the detachment of the vortices, the os- 
cillating sphere still sustains the vortex tangle. 

Figure[3]shows the vortex line length at different oscillation 



4 




time (ms) 

FIG. 3: (Color online) Vortex line length at different velocities of the 
oscillation. (1) 30 mm/s, (2) 50 mm/s, (3) 70 mm/s, (4) 90 mm/s, and 
(5) vortex line length in the absence of the sphere. 



velocities as a function of elapsed time. Only the vortex line 
length inside the computational box (40 /jm)^ is calculated. In 
Fig.|3] the line length increases in the first 0.5 ms as the vortex 
rings are injected. After the vortices are attached, the vortices 
are stretched and correspondingly the line length continues 
to increase. However, when detached vortices escape from 
the computational box, the line length suddenly drops. The 
loss of the vortices balances out the injection and the growth 
of the vortices, so that the line length saturates. Although a 
slight increase in the line length can be seen for a velocity 
magnitude of 30 mm/s, the line length almost lies on line (5) in 
Fig. [3] which means that the vortices are not stretched by the 
sphere. It can be seen that the saturated value of the line length 
increases as the oscillation velocity magnitude increases. For 
velocity magnitudes above 50 mm/s the saturated line length 
value is larger than the injection of vortices, which suggests 
that vortex tangles are forming around the sphere. 



IV. DRAG FORCE 

In this section, we propose a method to numerically evalu- 
ate the drag force exerted on a sphere oscillating in quantum 
turbulence, and calculate it as a function of the magnitude of 
the oscillation velocity. 

A sphere with acceleration dUp/dt experiences a force ex- 
erted by the fluid given by^ 



2 dUr, 



(6) 



Furthermore, in quantum turbulence, an extra drag force acts 
on the sphere for the following reason. In the equilibrium 
state of quantum turbulence, the injected vortex rings are con- 
stantly grown by the sphere. The sphere pushes the fluid and 
does work on the fluid because the kinetic energy of the fluid 
increases due to the stretching of the vortices, and a conse- 
quent reaction force acts on the sphere. The force described 
above is the order of 0.01 pN with the parameters used in this 
paper, being found negligible compared with the additional 
quantum drag force 1 pN, obtained later of the order of 1 pN. 



computational box 



V 



FIG. 4: (Color online) Schematic view of the computational box and 
volume V. 



By considering the energy given to the fluid, the drag force 
can be estimated. The derivative of the kinetic energy K of the 
fluid with respect to time is given by 



dK 



1 



dt dt Jv 2^^ ^ 



(7) 



where the volume V surrounds the sphere and is chosen to be 
large enough so that the superfluid velocity vanishes on the 
outer surface of the volume (Fig.|4]i.^^ Note that the derivative 
and the integral are not interchangeable because the domain of 
integration changes every time the sphere moves. To handle 
the difficulty, we extend the domain of integration to inside 
the sphere and make the domain time independent. Instead 
we define the density of the fluid as a function of time and 
position. 



p{r,t)^PsQ[\r~x{t)\^a], 



(8) 



where 9(x) is a step function, x{t) is the central position of 
the sphere, and a is the radius of the sphere. From Eq. ([8]l and 
Euler's equation, Eq. ^ becomes 



dK 
'dt 



Kn 2 dt 



dV 



Vail 



1 2^? „, 



pir,t)vs-^dV 

Vail 

dvs 



1 2^? 
Vail 2 dt 



- f Vs {- Vp - p,{vs ■ V)vs} dV 
Jv 



I -v,^^dV- [ p,Vs{vs-V)vsdV 
Jv^n 2 ot Jv 

- f Vs-VpdV, 

Jv 



(9) 



where the volume Van is the sum of V and the volume inside 
the sphere. 

We derive that the first and second terms on the right-hand 
side of Eq. (|9]) are canceled out. The time derivative of the 



5 



integrand in the first term becomes 
I = p.|e[|r-a.(f)|-«] 



p,6[|r-x(f)|-a](V^|r-a;(f)|) 



dx{t) 
dt 



= -p,5[|r-a:(0|-a]^-f|j|-M0, (10) 

where V;c is the gradient with respect to x. The delta function 
in Eq. ( fTOl l transforms the integral of the first term on the right- 
hand side of Eq. (|9]l into the surface integral, 

/ \v}%dV = - / ip,v,25(|r-a;|-fl) 
Jk,, 2 dt Jk,, 2 



all 
1 X 

\r — x\ 



UpdV 



f 1 2 



where the area S is the surface of the sphere, and dSx points 
outside the sphere. 

On the other hand, the incompressibility V • = leads 
to Vs{vs ■ 'V)vs — ^ijv/vs) and allows us to rewrite the 
second term on the right-hand side of Eq. (|9]) 

= [ psl;Vs^Vs-dS, 
Js' 2 

where S' is the sum of S and the outer surface of V . Since the 
velocity on the outer surface of V vanishes, only the integral 
over the sphere remains, and one obtains 

J^PsVs{vs-V)vsdV = - J^^PsV,^Up-dSi. (11) 

Thus the first and second terms on the right-hand side of Eq. 
^ are canceled out. Then Eq. ^ becomes 



dK 
dt 



= - / Vs- VpdV 
Jv 

= - [ V-{pVs)dV 
Jv 

= -j^pVs-dS2, (12) 



where the incompressibility is used, and dS-2 points inside 
the sphere. From Eq. (|4|i, Vg ■ dS-2 = Up ■ dS2, and Eq. (fTST i 

becomes 

^ = - J^pUp dS2 = -Up - (j^pdS2j ■ (13) 



In Eq. (fT3T l. the integral equals the force 
fluid on the sphere, and one obtains 



dK 
dt 



Up ■ Ff. 



exerted by the 



(14) 



The force Ff^s is resolved into the drag force .Fdrag parallel 
to Up and the lift force F^ft perpendicular to Up, and finally 
Eq. ( fT4b becomes 



dK 
dt 



Up ■ -fdrag- 



(15) 



Here we assume that Fdrag is antiparallel to Up averaged 
over a time period Tint much longer than the oscillation pe- 
riod, and that Up and i^drag oscillate in phase. The first as- 
sumption is based on the fact that a drag force giving energy 
to the surrounding fluid works against motion of an object. 
For the second assumption, it is considered that the drag force 
works when the vortices on the sphere are stretched, and the 
numerical simulation of Fig. |2] clearly shows that the vortex 
stretch is caused by the motion of the sphere. As a result, the 
drag force is synchronized with the oscillation of the sphere 
through the vortex growth without delay, and the second as- 
sumption is found to be valid. Then, taking the average of Eq. 
over Tint, we obtain 



dK 
dt 



Up ■ -fdr: 



■ag/ 



pQ 
p drag ' 



(16) 



where Wp and F^^^^^ are the amplitude of each oscillating quan- 
tity. Consequently the drag force is evaluated by dividing 
(dK/dt) by the velocity of the sphere. In the equilibrium 
state of quantum turbulence, the increase in the kinetic energy 
of the fluid in volume V equals to the energy of the vortices 
that escape from the computational box.^^ Thus (dK/dt) can 
be evaluated by averaging the energy of vortices that escape 
from the box over time period Tint. The energy of vortices 
can approximately be evaluated from its total line length. The 
energy of the vortex per unit length is given by 



4% 



In 



(17) 



where L is the characteristic length in the system, namely the 
computational box in this case. The line length of the vortices 
escaping from the box is measured over Tint = 100 ms, and the 
energy is calculated by multiplying Eq. ( fTTI l by the length.^— 
The line length also includes that of the injected vortices so 
the length of the injected vortices is subtracted from that of 
vortices that escape from the box. 

The drag force is plotted in Fig.|5]as a function of the os- 
cillation velocity. It should be noted that the unit of the drag 
force is pN, which is hundredths of the previously presented 
experimental data, say, 0.1 nN in Ref. 10 and 0.05 nN in Ref. 
14. This is explained by the difference in the object sizes used 
in the present study in comparison with these previous exper- 
imental studies. In Ref. 10, Jager et al. used an oscillating 
sphere with radius 100 /jm, which is a hundred times larger 
than ours, and accordingly the vortex tangle created around 
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10 100 40 60 80 100 

velocity (mm/s) velocity (mm/s) 

FIG. 5: The drag force and drag coefficient as a function of the ve- 
locity. The dotted line on the left-hand figure indicates the slope of 

the sphere would be larger than that of the present results. In 
Ref. 14, the length of wire over which vortex tangle is created 
is not known because the legs of the wire are installed on the 
experimental cell and only the middle of it can vibrate. Let 
us suppose a tenth of the wire effectively contributes to the 
turbulence. The diameter of the wire is 2.5 /jm and the whole 
length is 2 mm (then 200 /jm of the wire creates the turbu- 
lence), and the ratio of the wire size to the sphere is about 
60. It is expected that the wire could create at least 50 times 
of the size of the vortex tangle of our results. Hence, the ob- 
served difference in magnitudes of the drag force found by the 
present study and the previous experimental observations may 
be reasonable. 

In the range of velocity magnitudes from 50 to 90 mm/s, 
the drag force agrees well with a profile. Substituting the 
values of the drag force to Eq. ([T]i yields the drag coefficient, 
and such results are plotted in Fig. |5] The order of the drag 
coefficient in Eq. ([T]i is unity, being consistent with the values 
found in the related experiments. 

Skrbek and Vinen discussed in Ref. 17 that quantum turbu- 
lence created by oscillating objects at zero temperature also 
exhibits features similar to those of classical counterpart by 
citing various experiments. As we mentioned in Sec. HI the 
similarity between these phenomena is confirmed in homoge- 
neous and isotropic turbulence, as the form of Kolmogorov's 
law. It is remarkable that this similarity is numerically con- 
firmed in another type of turbulence, say, turbulence created 
by oscillating objects. 

Skrbek and Vinen'^ also described the existence of two pos- 
sibilities for the flow of the superfluid. One is quasiclassical 
laminar flow, in which a vortex tangle of low density is formed 
around the objects and the flow of the superfluid mimics a 
classical laminar flow. The other is quasiclassical turbulent 
flow, in which the density of vortex tangle becomes high and 
large scale rotational flow appears, mimicking classical tur- 
bulence. From Fig. |2] vortices of large scale comparable to 
the sphere size can be seen, and in fact, the flow exhibits the 
behavior of Eq. ([T]i. These agreements evince the similarity 
between the classical and quantum turbulences and support 
the picture of quasiclassical turbulent flow. 

The process whereby the turbulence is created raises the 



TABLE I: The drag coefficient computed for different time intervals. 
Oscillation velocity is 50 mm/s. 



T(ms) 


0.03 


0.05 


0.1 


Cd 


0.78 


0.68 


0.59 



possibility that the time interval i in which the vortex rings 
are injected can affect the drag coefficient. Table|T]shows the 
drag coefficient calculated using different time intervals at an 
oscillation velocity of magnitude 50 mm/s. As the injection 
interval becomes shorter, the rate of vortex growth increases, 
resulting in an increase in the drag coefficient. However, it 
may be noted that the drag coefficient remains of order unity. 

V. CONCLUSION AND DISCUSSION 

Quantum turbulence was created around an oscillating 
sphere, and its evolution was studied. We proposed a method 
to evaluate the drag force in quantum turbulence. The drag 
force was calculated by considering the energy supplied by 
the sphere, and it is proportional to the square of the magni- 
tude of the oscillation velocity. This dependency on the ve- 
locity is quite similar to that in classical turbulence. The drag 
coefficient was also calculated, and its order of magnitude is 
in agreement with that of the previous experimental results. 
The similarity between classical and quantum turbulences was 
confirmed by numerical simulations. 

We believe that the drag coefficient has an universal value 
independent of the details of the system. Our future work is 
to confirm whether the drag coefficient still remains of the 
order of unity even if the parameters such as the frequency of 
the oscillation and the geometry of the oscillating object are 
changed. 

The simulations were performed for a smooth sphere, but 
the objects used in the experiments have surface roughness. 
The roughness may increase the drag force and the drag coef- 
ficient; when a vortex is detached from the sphere due to the 
oscillation, it is likely to leave a small vortex bridge over the 
pinning sites on the surface. The small vortex plays a role of 
"seed" for the growth of the vortex, and this is equivalent to 
injecting more vortices to the sphere, resulting in an increase 
in the drag force and the drag coefficient. However, there is 
no clear way to proceed to obtain the solution of the equations 
of motion of vortices under the boundary conditions which 
describe the surface roughness at present, and the numerical 
simulations require a model that phenomenologically depicts 
the vortex pinning. 

Figure |5] enables us to discuss the transition from the tur- 
bulent to laminar state. Below oscillation velocities of mag- 
nitude 50 mm/s, the drag force starts to deviate from Eq. ([1]) 
and tends to zero. This is because the motion of the sphere is 
so slow that vortices cannot be stretched before the attached 
vortices leave the sphere, resulting in the drag force tend- 
ing to zero. This means that the oscillating sphere can no 
longer sustain turbulence, and that the flow regime will re- 
turn to a laminar state. From dimensional analysis, the crit- 
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ical velocity at which turbulence returns to laminar flow is 
given by — cy^KoS, where c is constant with the order of 
unity, and this formulation agrees with the critical velocity in 
the experiments.^' In the simulation, K = 10^^ m^/s, co = 10^ 
rad/s, and finally one obtains vv ^ 30 mm/s. This is very close 
to the velocity 50 mm/s below which the deviation from Eq. 
([U starts in the simulation. Therefore the velocity magnitude 
50 mm/s can be considered in the simulations to represent the 
critical velocity magnitude. 
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